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ABSTRACT 

A comparison  via  Monte-Carlo  simulations  is  made  of 
Sliding  Mode  Observers  versus  Kalman  Filter  in  the 
homing  missile  guidance  system  using  different 
guidance  laws.  It’s  shown  that  the  sliding  mode 
observer  contributes  to  less  miss  distance  due  to  noise 
than  Kalman  filter.  Application  of  sliding  mode 
estimators  in  the  homing  loop  with  Proportional 
Navigation  guidance  and  with  a phase-lead 
compensation  based  on  sliding  mode  estimators  for  the 
flight  control  system  phase-lag  is  demonstrated.  It  is 
shown  that  this  system  gives  performance  comparable 
to  the  optimal  guidance  law  with  Kalman  filter  that  has 
perfect  knowledge  of  all  additional  information  it 
requires  for  estimation. 

INTRODUCTION 

An  ultimate  performance  criterion  of  a homing  missile 
guidance  system,  the  miss  distance,  is  crucially 
dependent  on  each  element  of  guidance,  navigation,  and 
control  of  a homing  interceptor.  In  particular, 
performance  of  the  estimator  (filter  or  observer)  of  data 
necessary  to  compute  a guidance  command  is  an 
important  factor  to  affect  the  miss  distance.  This 
performance  is  difficult  to  analyze  under  the  influence 
of  the  following  factors:  different  effects  caused  by 
parasitic  intercoupling  of  navigation  and  flight  control 
systems,  the  signal  phase  lag  due  to  bandwidth  limits  in 
navigation  and  flight  control  systems,  guidance 
command  saturation  in  presence  of  evasive  target 
maneuvers,  and  higher  measurement  noise  due  to  other 
target  countermeasures.  Estimator  contributes  to  each 
of  the  aforementioned  effects  in  a peculiar  way.  In  this 
case  numerical  Monte-Carlo  simulations  using  even  a 


simplified  mathematical  model  of  the  homing  loop  may 
give  a true  relative  picture  of  miss  distance  projections 
to  different  engagement  conditions  under  different 
estimation  procedures. 

In  this  paper  numerical  Monte-Carlo  simulations  are 
selected  as  a primary  instrument  to  compare  the 
performance  of  the  widely  used  Kalman  filters  in  the 
homing  loop,  presented  by  linearized  state-model  for 
the  relative  separation  between  missile  and  target 
perpendicular  to  the  fixed  reference1,  versus 
performance  of  Sliding  Mode  observers  based  on  an 
emerging  technique  in  the  Sliding  Mode  Control  (SMC) 
theory  - higher  order  sliding  modes2'3.  A unique  feature 
of  these  new  observers  is  that  they  separate  signal  from 
noise  (i.e.,  reconstruct  a signal  and  its  derivative  from 
noisy  measurements)  on  the  basis  of  a specified  limit 
for  the  signal  second  time-derivative. 

The  main  goal  of  this  experiment  is  to  study  robustness 
to  noise  (miss  distance  standard  deviation)  and  to  the 
target  step-constant  maneuver  at  different  times-to-go 
(mean  miss  distance)  of  the  guidance  system  with 
Kalman  filter  (KF)  and  with  SMC  filter/differentiator. 
These  two  criterions  ought  to  show  how  good  two 
estimators  track  given  signal  and  attenuate  noise 
relatively  to  each  other.  Separation  principle  for 
estimation  process  in  this  nonlinear,  overly  constrained 
problem  is  not  valid,  so  evaluating  performance  of  two 
algorithms  in  terms  of  estimation  error  without 
considering  homing  loop  does  not  make  a good 
projection. 

The  main  expectation  of  the  intended  experiment  is 
that,  having  less  input  information  (line-of-site  (LOS) 
only),  a SMC  robust-to-noise  filter/differentiator  will 
demonstrate  the  overall  guidance  system  performance 
comparable  to  that  with  KF,  which  require  in  addition 


n This  work  was  supported  by  the  US  Army  SMDC,  under  contract  number  DASG60-94-C-0045 
Distribution  Statement  A:  Approved  for  public  release;  distribution  is  unlimited. 


3-1 

UNCLASSIFIED 


UNCLASSIFIED 


to  LOS  measurements  the  following  data:  range  to 
target,  time-to-go,  missile  acceleration,  maximum 
expected  target  acceleration,  and  noise  standard 
deviation.  In  simulations,  we  assume  noise  in  LOS 
measurements  only  and  perfect  information  of  other 
parameters  necessary  for  KF,  as  well  as  perfect 
information  of  closing  velocity.  In  this  case,  KF 
performance  in  the  homing  loop  will  be  more  optimistic 
than  that  of  SMC  filter  as  far  as  real  world  situation  is 
concerned.  Thus,  demonstrating  comparable 
performance  in  this  idealized  test  will  indicate  clear 
advantage  for  SMC  algorithm. 

The  main  reason  to  carry  out  this  experiment  is  in  the 
following.  Quite  a few  works  has  appeared  recently  in 
the  literature4  8 studying  application  of  SMC  methods  to 
homing  missile  guidance  and  flight  control  systems 
design.  Theoretical  and  numerical  results  look 
promising.  There  were  shown  clear  advantages  in 
robustness  and  the  missile/target  acceleration  ratio.  The 
SMC-based  navigation  data  processing  can  be  made  the 
last  step  towards  designing  a truly  integrated  homing 
missile  guidance,  navigation,  and  control  system  within 
a unified  framework  of  SMC  design  methods  including 
recent  advances  in  the  field  of  higher  order  sliding 
modes.  Thus,  new  principles  in  estimation  and  control 
decisions  for  missile  guidance,  navigation,  and  control 
are  really  worth  studying  and  evaluating. 

PRELIMINARIES 


Guidance  System  Model 

As  known1,  the  linearized  guidance  system  model 
serves  as  a good  approximation  of  nonlinear 
engagement  kinematics  for  the  purpose  of  comparative 
analysis  of  different  guidance  laws.  True  performance 
projections  are  obtained  even  under  nonlinear 
constraints  (such  as  acceleration  command  and 
command  rate  saturation)  if  the  engagement  happens  in 
the  vicinity  of  the  collision  course  and  the  closing 
velocity  is  relatively  constant  during  the  flight 
(excluding  the  last  moment,  when  the  rapid  reverse  in 
the  sign  of  the  closing  velocity  occurs). 

The  structure  of  a linearized  guidance  system 
presented  in  the  work1  will  be  used  to  investigate  the 
performance  of  the  system  in  presence  of  measurement 
noise,  and  using  different  estimates  of  Line-of-Sight 
(LOS)  rate.  The  following  constraints  on  missile 
acceleration  and  acceleration  jerk  will  be  used 
|rit.|<30G,  |«^|  < 100G/S. , and  implemented  as  the 

“Performance  Limits”  block  given  in  Fig.  1 where  the 
flight  control  system  is  modeled  as  a single-lag  block 
with  the  time-constant  T =0.5 S. 


Fig.  1 “Performance  Limits”  Block 
Then,  the  guidance  system  will  have  the  structure 


Fig.  2 Guidance  System 

The  linearized  state  equation  for  the  relative 
separation  between  missile  and  target  perpendicular  to 
the  fixed  reference  is 

y = nT-nL  (1) 

where  the  target  acceleration,  rij , is  modeled  as  a 
constant  maneuver  with  6G  amplitude,  such  that 
acceleration  ratio  is  relatively  large  but  finite, 

Hc  UM  _5 

H'r 
1 max 

The  relative  separation  along  the  fixed  reference  is 
supposed  to  obey 

x = Vc{tF~t)  (2) 

where  the  closing  velocity  Vc  =9000  Ft/S.  is  held 
constant,  and  the  homing  time  tF  =10S.  if  not  posted 
otherwise.  It’s  assumed  that  the  missile  initially  is  on  a 
collision  course,  and  that  the  target  evasive  step- 
constant  maneuver  will  occur  at  different  times-to-go 
before  the  collision.  The  expected  peak  miss  distance 
will  be  at  the  time  which  is  large  enough  to  escape 
safely  but  too  small  for  the  tracking  filter  to 
accommodate  a sudden  change  in  the  signal  behavior. 

The  measurement,  LOS  angle  A,  is  sensed  with 
noise 

As  =A  + An  (3) 

where  An  is  a Gaussian  white  noise  with  zero  mean  and 
standard  deviation  of  10~^  Rad  (l  mR). 

Line-of-sight  rate  (LOS),  CO % , is  estimated  from 
measurements  As  by  the  three-state  discrete  Kalman 
filter  presented  on  p.166  in  the  reference1  and  by  a 
SMC-based  robust-to-noise  filter/differentiator 
presented  in  Section  2.3. 
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Different  guidance  laws  will  be  simulated  starting  with 
Proportional  Navigation  law, 

nc  = N'Vcd , CO x = A . (4) 

The  expectations  are  that,  having  less  information 
requirements,  SMC  observers  can  keep  good  balance 
between  phase  lag  in  the  LOS  rate  estimate  and  noise 
attenuation,  and  achieve  performance  comparable  to  an 
idealized  KF  (with  perfect  information  about  the  range- 
to-target,  Rjm  ; time-to-go,  tgo  ; missile  acceleration, 

n i ; maximum  expected  target  acceleration  iij  max  , 
and  noise  standard  deviation.) 


Kalman  Filter  in  the  Homing  Loop 


The  tree-state  (y,  y,  hj  j digital  Kalman  filter  i 
used  to  estimate  LOS  rate  6) A , 

y + yfgo 


IS 


coA  =■ 


V t2 

v cl  go 


and  target  acceleration,  hj , which  is  identical  to  the 
one  shown  in  Fig.  9.2  in  the  reference1,  except  that  the 
achieved  missile  acceleration,  , rather  than  the 


commanded  acceleration,  nc  , is  fed  back  into  the  filter. 
Different  modifications  of  the  program  code 
“C9L3.cpp”,  developed  in  the  reference1.  Chapter  9,  are 
used  to  simulate  the  guidance  system  in  Fig.  3 with  KF 
and  different  guidance  laws. 


SMC  Filter/Differentiator  in  the  Homing  Loop 

Developed  recently3,  an  exact  differentiator,  based  on 
the  second  order  SMC  method,  can  reconstruct  a given 

signal  fit ) and  its  derivative  fit ) provided  certain 


: J"  x 2<ir, 


*2  = Pi 


(f(t)-x  i) 


+ Po  J sgn(/ (0  _ X\)dr. 


(5) 


0.5 


If  pQ  >4L,  pl  >0.5 L 

then  Xj  — > fit),  xA  — > f(t)  in  a finite  time,  as 

proved  in  the  work3.  Since  the  structure  (5)  uses  the 
signal  state-model 

*1  =*2- 


X2  = fU X \f(t)\ 


<L, 


to  estimate  the  states  Xj  ,x2  of  a signal,  it  can  be  called 
the  state-observer. 


In  case  of  the  violation  of  condition  |/(0|  — L , say, 

fit)  = f0  it)  + vit),  | f0  (f)|  < L,  |v(0|  » L,  (6) 

the  observer  (5)  will  reconstruct  the  actual  signal 
fQ  it)  with  good  accuracy,  Xj  — > f0  it) ; however, 
x2  it)  will  have  significant  high-frequency  component 
due  to  noise,  v(t) , amplification  on  the  first  “high- 
gain”  term  in  the  expression  for  x2(0  (Eq.  (5)). 
Moreover,  because  x2  it)  cannot  follow  too  aggressive 
fit) , in  digital  implementations  of  the  observer  (5), 
there  will  appear  so-called  “chattering”  of  term  x2  (t) 
due  to  a finite  time-step  of  the  algorithm.  However, 
Xj  (/)  will  be  smoothening  out  by  an  integrator.  To 

estimate  f0  it) , one  can  feed  x2  (/) , which  contains 
information  about  / 0 it) , to  the  input  of  another 


Fig. 3 Robust-to-Noise  Differentiator 


constraint  on  the  signal  second  derivative  /(f)  <L, 
for  known  L > 0 . It  has  the  following  structure 


observer  of  the  form  (5),  and  reconstruct  f 0 it) , 
applying  certain  constraint  Lj , fQ  (f )J  < Lj . The  value 
of  Lj  should  be  a result  of  a tradeoff  between  the  ripple 
magnitude  and  the  phase  lag  in  f 0 ( t ) estimate. 
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Thus,  for  noisy  signals  satisfying  (6),  two  observers 
connected  in  series  (such  that  the  output  x2  (/)  of  the 
first  one  is  fed  into  the  input  of  another,  and  the  state 
Xj  (?)  of  the  second  one  is  the  designated  output,  as 
shown  in  Fig.  3,  become  a robust-to-noise 
differentiator.  This  system  will  separate  a desired  signal 
from  noise  by  the  factor  of  the  actual  signal  second 
derivative.  This  important  time-domain  property  of  a 
signal  makes  this  observer/differentiator  to  be  unique 
among  various  filters  and  state-observers,  which  use 
statistical  or  spectral  signal  properties. 


The  second  order  sliding  mode  scheme  presented  in  (5) 
is  known  as  the  Super  Twisting  Algorithm  developed 
by  Levant3.  Another  SMC  estimator,  based  on  Twisting 
Algorithm  by  Yemelyanov,  Korovin,  Levantovsky2'3,  is 
working  as  a filter  producing  estimates  of  a given  signal, 
x(t) , and  its  derivative,  x(t) , using  their  noisy 
measurements 

Tl(0  = x(0  + V1(0, 

y2(t)  = x(t)  + v2(t). 

A robust  2-sliding  estimator  based  on  Twisting 
Algorithm  is  given  as 

x = Pi  sgn(y2  ~ x)  + p0  sgn(^1  - x) , (8) 


In  absence  of  noise,  2-sliding  mode  in  the  system  (8) 
provides  for  x — > x(t),  X — > x(t)  in  a finite  time  if 
|x(f)|<L.  Given  noisy  signals  y\(t\y2(t) , the 
high  frequency  switching  of  Sgn(e)  and  Sgn(e)  will 
reconstruct  x(t ) and  x(t)  to  the  extent  where 

\x(t)\  < L . So,  if  \x(t)\  < L and  Iv'j^  (0 1 » L . 


Then  we  have  a good  x(t),x(t)  recovery  in  a 
practical  filter  of  the  form  (8). 


The  structure  of  the  LOS  rate  estimator  based  on  the 
observer  (8)  is  presented  in  Fig.  4. 


The  last  structure  will  be  used  for  LOS  rate  estimation, 
and  the  Super  Twisting  Algorithm  will  be  employed 
further  for  the  purpose  of  compensation  for  flight 


control  system  dynamics  via  a phase-lead  block  based 
on  the  structure  in  Fig.  3. 

Thus,  the  following  equations  of  Twisting  SMC 
estimator  together  with  guidance  algorithm  and  flight 
control  system  dynamics  are  used  in  simulations 


Z\  =60(z2  -Zi), 
z2  = 60(Ay  - z2 ), 


Analog  Prefilter, 


A[A:-7}]=500zj[&-7^],  SampledLOS, 


Z3  - Z4 , 

Z4  = 25  • (sgn(  X[k  ■ Ts  ] - z3)  + 


sgn( 


X[k  -Ts]-X[(k  -D-T,] 


- ^4 


SMC-differentiator 

r .1 


z5  =50|  3|z4  -z5|2  -sgn(z4  -z5)  + 4z6 

v 

z6  =sgn(z4  -z5), 

SMC-filter 

C02  = 0.002  ■ Z5 , LOS  rate  estimate, 
nc  = N'Vcd)2 , N ’=  3;  4, 

Proportional  Navigation  Guidance  Law, 
nc\k]  = SA T966 (nc [k  • Ts  ]), 

= ~SA ?3220  (2>  (nL  ~nc[k ])) 

Performance  Limits  Block  (T  = 0.5 sec ) 
L ■ sgn(x),  |x|  > L, 
x,  IjcI  < L. 


Function  SAT^{x)  = 


966  -^r  equivalent  to  30G, 


sec 


3220  — equivalent  to  100  G/S. 
sec3 


(9) 
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A digital  computer  in  the  homing  loop  processing 
measurement  information  and  calculating  the  guidance 
command  is  simulated  as  given  in  Fig. 5. 


flight  times  comparable  to  the  effective  guidance 
system  time-constant.  Optimal  guidance  law  to 
compensate  for  a single  phase  lag  in  flight  control 
system  has  been  derived1, 


Fig. 5 Discrete  Part  of  Guidance  System 


The  sensed  continuous  signal  As  being  low-pass  pre- 
filtered is  coming  to  the  A/D  converter  represented  by 
S&H  (Sample  & Hold)  element  with  the  sampling  time 
ts  =0.01Sec.  Then,  the  stepwise  constant  digital 

signal  A\k-Ts]  is  used  by  a discrete  guidance 
algorithm  to  calculate  nc . The  guidance  command  is 
stepwise  constant,  nc  [k  ■ Ts  ] , and  it  is  updated  with 
the  same  rate  on  another  S&H  element  representing 
D/A  converter.  When  a discrete  Kalman  filter  is  used  to 
produce  necessary  estimates,  the  step  of  the  digital 
computer  algorithm  for  a discrete  Kalman  Filter  is 
equal  to  the  sampling  time.  In  case  of  SMC-based 
differentiators  and  filters  calulating  COA , the  rate  of  the 
digital  computer  algorithm  for  discrete  implementation 
of  SMC-estimators  is  higher  than  the  sampling  rate 
(SMC-estimators  are  dynamic  between  samples). 

Advanced  Guidance  Laws  based  on  Kalman  Filter 
Estimation  and  Proportional  Navigation  Guidance 
with  Phase-Lead  Network  based  on  SMC 
Differentiators 


One  significant  difference  between  Kalman  filter  and 
SMC  differentiator  is  that  Kalman  filter  is  a three-state 
estimator,  while  SMC  filter/differentiator  estimates 
only  two  states:  LOS  and  LOS  rate.  The  third  estimates 
in  Kalman  filter,  target  acceleration,  is  used  to 
implement  more  advanced  guidance  laws.  Therefore, 
fair  comparison  should  be  made  of  any  available 
guidance  law  using  Kalman  filter  with  any  guidance 
law  that  can  be  implemented  using  SMC  differentiator 
only.  It’s  known1  that  Proportional  Navigation  (PN)  is 
optimal  for  zero-lag  guidance  system  in  absence  of 
target  maneuvers,  while  Augmented  Proportional 
Navigation  (APN)  is  optimal  for  zero-lag  guidance 
system  in  presence  of  step-constant  target  maneuver 

nc  =N,fycC0A  + 2,1t)' 

However,  phase  lag  in  flight  control  system  and  the 
estimator  itself  is  a source  of  significant  miss  at  short 


nc=N  1 


N'= 


VC(0A  + 0.5nT  -nL  ■ 

xl 


(e~x  +x-l) 


6x2  ( e x + x — 1) 


2x3  + 3 + 6x  - 6x2  - \2xe  x — 3e" 


lx 


, X- 


1 go 
T ' 
(11) 


Thus,  having  additional  information  about  n tg0 
and  using  estimates  COA , fij  from  Kalman  filter,  one 
can  implement  guidance  laws  (10)  and  (11). 

The  main  concern  of  the  optimal  guidance  law  (OGL) 
(11)  is  to  compensate  for  parasitic  single  lag  dynamics 

nL  = — ~ — nc  • Knowing  about  a phase  lag  in  flight 
1 + sT 

control  system,  one  can  design  a lead  block,  1 + sT , to 
compensate  for  known  phase  lag,  using  SMC  robust-to- 
noise  differentiator  presented  in  Fig.  3.  Using  the  lead 
block  1 + sT  , we  can  expect  significant  reduction  of 
mean  miss  at  short  flight  times,  although,  for  the 
expense  of  increase  in  miss  standard  deviation  since 
differentiation  leads  to  amplification  of  noisy 
component.  So,  an  alternative  to  OGL  (11)  may  be  the 
PN  guidance  law  with  a phase  lead  based  on  two 
cascades  of  SMC  differentiators.  The  idea  is  to  use 
SMC-differentiator  one  more  time  to  produce  LOS 
second  derivative  estimate  and  create  LOS  rate  estimate 

with  a phase  lead  as  &A  ieacj  = (bA  + 0.5^  . Then, 
^Alead  t>e  used  to  compute  the  PN-guidance 
command. 


Thus,  PN  guidance  law  with  a phase  lead  is 
implemented  as  follows.  Equations  (9)  are  amended 
with 


(10) 
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£7  ~;70’ 

z8  =sgn(z5  -z7), 

\no={' 


8|z5  -z7 1 2 -sgn(z5  -z7)  + 15z8 


V 


z9  -3 


%o  - Z91 2 -sgn(^0  -z9)  + 15zio 


Z10  = sgn(^0  - z9 ), 


= 0.002  ■ Z5 , 

= 0.002  • z9 , 

lead  = $>A  + °-5  ' ^2  > 
nc  — N VcCOX  lead’  N — 4. 


(12) 

Thus,  another  important  round  of  numerical  tests  is  to 
compare  performance  of  APN  (10)  and  OGL  (11)  based 
on  Kalman  filter  estimates  with  the  work  of  PN  with  a 
phase  lead  based  on  two  cascades  of  SMC  robust-to- 
noise  differentiators  (9),  (12). 


EXPERIMENT  OBJECTIVE  AND  SETUP 

A flight  control  system  phase  lag  plus  a filter  phase  lag, 
plus  a phase  lag  introduced  by  a sampler  altogether 
with  the  commanded  acceleration  limit 
| nc  | max  = 30g  produce  a peculiar  miss  distribution 

pattern  under  the  presence  of  target  maneuvers  and 
measurement  noise.  We  believe  that  the  homing  loop 
model  and  the  selected  constraints  will  highlight 
performance  of  both  estimators  the  way  we  are  looking 
for. 

The  main  objective  of  the  current  experiment  is  to 
obtain  via  50-runs  Monte-Carlo  statistics  for  each 
particular  flight  time  within  range  [0.5 ;10  .OjSec  , 
for  a total  of  1000  times,  mean  miss  distance  and  miss 
distance  standard  deviation,  using  with  both  estimators 
and  different  guidance  laws  discussed  in  Section  2. 


flight  times  due  to  overall  time  lag  in  the  guidance 
system  and  actual  instability  of  the  homing  loop. 

Computer  simulations  were  accomplished  ug^o^,^  gjyjQ 
Microsoft  Visual  C++  5.0  and  the  program  code 
“C9L3.cpp”  developed  by  P.  Zarchan1  for  a PN 
guidance  system  with  Kalman  Filter.  There  were  made 
modifications  of  this  program  code  for  the  guidance 
system  with  SMC  estimator  amd  different  guidance 
laws. 

SMC -filter 

A case  was  considered  in  which  there  was  a constant  6- 

g target  maneuver  and  1 miliradian  of  measurement 

noise  <7 nolse  = 0.001/tarf  . Initial  Heading  Error  is 

zero.  Kalman  filter  has  absolute  knowledge  of 

ma  x(«r  >’  ho  • RTM  > nL  and  noise  stanfjard^^,  Nayjg 

deviation.  The  sampling  time  of  measurement  and 

command  update  is  Ts  = 0.01.5  (Kalman  filter  works 

at  this  sampling  rate).  The  integration  time  step  is 

T = 0.000 IS  (the  continuous  part  of  the  homing  loop 

and  the  SMC  filter  work  at  this  sampling  rate).  The 

program  “C9L3.cpp”  uses  Euler  integration  method. 

MONTE-CARLO  SIMULATIONS 


First,  we  simulate  PN  guidance  law  with  the  effective 
navigation  gains  N’=  3 and  N’=  4 and  with  LOS  rate 
estimate  by  Kalman  filter  and  by  SMC  differentiator. 
The  goal  is  to  observe  the  sensitivity  of  the  guidance 
system  to  a step  constant  target  maneuver  with 

acceleration  ratio  = 5 and  measurement  noise 

1lT 

'max 

in  the  presence  of  a single-lag  flight  control  system. 
The  results  of  simulations  are  given  in  Figs.  6-9. 

The  next  set  of  simulations  is  to  identify  what  possible 
advantage  SMC  differentiatiors  can  contribute  to 
compensation  of  a single-lag  flight  control  system  as 
compared  to  advanced  guidance  laws  using  Kalman 
filter  estimates.  Performances  of  APN  (10)  and  OGL 
(11)  based  on  Kalman  filter  estimates  and  PN  with  a 
phase  lead  based  on  two  cascades  of  SMC  robust-to- 
noise  differentiators  (9),  (12)  are  compared.  The  results 
of  simulations  are  given  in  Figs.  10,11. 


Certainly,  we  can  consider  that  the  mean  miss  is  mostly 
due  to  the  target  maneuver,  while  the  standard  deviation 
of  the  miss  is  due  to  noise1.  As  far  as  the  initial  heading 
error  is  zero,  we  can  interpret  the  flight  time  tp  as 
time-to-go  at  which  a constant  6-g  target  maneuver 
occur1.  We  can  expect  that  a large  miss  happens  at  short 


Short  summary  of  important  benchmarks  is  presented  in 
Table  1. 
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Table  1 


Guidance 

Law 

Estimator 

Peak  Mean 
Miss  (Ft) 

Steady -state 
Mean  Miss  (Ft) 

Peak  Miss 
Std.  Dev.  (Ft) 

Steady-state 
Std.  Dev.  (Ft) 

PN,  N’=3 

KF 

51 

0.3 

11 

1.0 

PN,  N’=3 

SMC 

54 

0.012 

0.7 

0.1 

PN,  N’=4 

KF 

34.5 

0.3 

10.8 

1.16 

PN,  N’=4 

SMC 

32.7 

0.01 

0.82 

0.11 

APN,  N’=4 

KF 

29.3 

0.25 

1.98 

OGL 

KF 

16.6 

0.06 

7.7 

1.05 

PN,  N’=4 
with 

Phase  Lead 

SMC  (9) 
and 

SMC  (12) 

14.4 

0.001 

2.5 

0.8 

Summary 

The  guidance  system  with  SMC  estimator  is  more 
robust  to  noise  than  with  the  optimally  tuned  Kalman 
filter.  However,  it  is  known  that  the  detuning  of 
Kalman  filter  may  improve  performance.  On  the  other 
hand,  performance  of  Kalman  filter  will  degrade  from 
the  values  obtained  in  this  experiment  if  there  is  noise 
in  the  additional  information  (range  to  target,  time-to- 
go,  missile  acceleration).  Particular  benefits  in  using 
SMC  methods  for  missile  guidance  are 

1.  The  SMC  estimators  are  proved  to  be  effective  as 
robust-to-noise  differentiators, 

2.  SMC  estimators  require  less  input  information  than 
KF  to  produce  a necessary  estimate. 

Application  of  SMC  estimators  in  the  homing  loop  with 
PN  guidance  and  with  a phase-lead  compensation  for 
the  flight  control  system  phase  lag  gives  performance 
comparable  and  actually  very  close  to  the  optimal 
guidance  law  with  Kalman  filter  and  perfect  knowledge 
of  all  additional  information  it  requires  for  estimation. 

CONCLUSIONS 

In  this  work  we  have  compared  the  performance  of 
Proportional  Navigation,  Augmented  Proportional 
Navigation,  and  Optimal  Guidance  Law  with  a three- 
state  Kalman  filter  versus  Proportional  Navigation  with 
and  without  a phase-lead  compensator  employing 
Higher  Order  Sliding  Mode  Observers.  It  was  shown 
that  one  could  expect  the  good  work  of  SMC  estimators 
in  a missile  guidance  system.  Given  less  information 
requirements  for  the  SMC-based  estimators  than  for 
Kalman  filters,  one  can  expect  benefits  in  applying  this 
technology  to  a homing  interceptor.  The  next  step  in 
this  research  effort  is  to  create  an  integrated  GNC 
system  for  an  advanced  interceptor  using  SMC  theory, 
and  derive  the  advanced  guidance  laws  which  guarantee 
robust  to  target  maneuvers  intercept  under  less 
acceleration  ratio  requirements. 


Flight  Time  (S) 


Fig. 6 Mean  Miss  versus  Flight  Time  for  PN  guidance 
with  KF  and  SMC  observer 
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Flight  Time  (S) 


Fig.7  Miss  Standard  Deviation  versus  Flight  Time  for 
PN  guidance  with  KF  and  SMC 


Fig. 8 Mean  Miss  versus  Flight  Time  for  PN  guidance 
with  KF  and  SMC  observer 


Flight  Time  (S) 


Fig.9  Miss  Standard  Deviation  versus  Flight  Time  for 
PN  guidance  with  KF  and  SMC 
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Fig.  10  Mean  Miss  versus  Flight  Time  for  APN,  OGL 
with  KF,  and  PN  with  SMC  observer  and  a phase-lead 
compensator 


Fig.l  1 Miss  Standard  Deviation  versus  Flight  Time  for 
APN,  OGL  with  KF,  and  PN  with  SMC  observer  and  a 
phase-lead  compensator 
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